Non-adiabatic elimination of auxiliary modes in continuous quantum measurements 
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When measuring a complex quantum system, we are often interested in only a few degrees of 
freedom — the plant, while the rest of them are collected as auxiliary modes — the bath. The bath can 
have finite memory (non-Markovian), and simply ignoring its dynamics, i.e., adiabatically eliminat- 
ing it, will prevent us from predicting the true quantum behavior of the plant. We generalize the 
technique introduced by Strunz et. al. [Phys. Rev. Lett 82, 180f (1999)], and develop a formalism 
that allows us to eliminate the bath non-adiabatically in continuous quantum measurements, and 
obtain a non-Markovian stochastic master equation for the plant which we focus on. We apply this 
formalism to three interesting examples relevant to current experiments. 



Introduction. — Recent developments in techniques of 
high-precision metrology have allowed quantum-level 
measurement and control of matters of at all scales, rang- 
ing from single atoms [ij to macroscopic mechanical oscil- 
lators [H. In these experiments, the atoms or mechanical 
oscillators, as objects of interest (or the plant), are usual 
coupled to auxiliary degrees of freedom (or the bath), e.g. 
the cavity mode in cavity QED systems, which are in turn 
coupled to external readout devices. It is often desirable 
to obtain a self-contained equation for the state of the 
plant, by eliminating bath degrees of freedom, especially 
when we want to implement a real-time feedback control. 
In the literature, the simplest approach is to ignore the 
dynamics of bath modes by assuming that they follow 
the plant instantaneously, and can be adiabatically elim- 
inated. However, this becomes inadequate when bath 
modes evolve at scales longer than the plant, i.e., when 
the system becomes non-Markovnian. 

One way to account for a non-Markovian bath is the 
Fcynman- Vernon influence functional method Q . Diosi 
and Strunz et. al. 0-0 developed an equivalent (but 
much simpler) method by unraveling the bath evolu- 
tion into possible quantum trajectories. These trajec- 
tories arc shown to drive a non-Markovian stochastic 
Schrodinger equation (SSE), which average into the exact 
non-Markovian master equation. Although their model 
does not include measurement a priori, the SSE at the 
Markovian limit can be interpreted as the evolution of 
the plant's pure state under continuous measurement. 
In general, however, the physical interpretation of non- 
Markovian SSE in terms of measurement has yet to be 
clarified, as discussed by Diosi and Wiseman et al. Q. 

Here we consider the non-Markovian measurement pro- 
cess involving a plant-bath system, in which the bath 
has finite memory, and measurement is done through 
the bath. We further assume that the incoming probe 
field for measurement is a quantum Wiener process ISl, 
and that the output field is projectively measured [10|. 
The setup is shown in Fig.[TJ In addition to the usual 
applied assumptions for the plant-bath interaction — the 
bath is bosonic and couples linearly to the plant, we as- 
sume that the bath is also linearly coupled to the probe 



field. We do not make assumptions on the plant, nor 
the plant quantity that couples to the bath. By general- 
izing the Diosi-Strunz approach, we show that the bath 
can be eliminated from the full evolution equation, re- 
sulting in a non-Markovian Stochastic Master Equation 
(SME) which governs the density matrix of the plant, 
and thereby has a distinctive physical meaning. 

The prescription we develop here can be applied to 
a wide range of non-Markovian quantum measurements, 
thereby laying the foundation for an interesting research 
direction. The purpose of this Letter is to present the 
general formalism and highlight three examples of rele- 
vance to current experiments, for which analytical forms 
of the SME can be obtained. Interestingly, in two of 
the three cases, the non-Markovian dynamics we obtain 
differs from conventional wisdom. 

Model. — The Hamiltonian for our model reads: 

H^Hp + Hh + Hint + Yj^^Vlkiak Hnit) + al knit)], 
Hb = ^JiUkalcLk, Hint = ^^hgk{Lal + Vau). (1) 

Here Hp, H^ and i^int arc the plant, bath, and interac- 
tion Hamiltonians, respectively; Ofe and Wfc are the anni- 
hilation operators and eigenfrequencies of different bath 
modes and [ofc, a\i] = Skk'', the plant operator couples 
to the bath through L, with gk its coupling constant to 
the A:-mode; bin{t) are annihilation operators for the in- 
put probe field at different times and [bin{t), ^in(^')] ~ 
S{t ~ t'); "fk is the coupling strength between the bath 
and the probe field. We exclude those modes that are 
not coupled to the probe field, as they will simply in- 




FIG. 1: (color online) Schematics of our measurement process. 
The plant is coupled to the bath, which in turn couples to 
an external probe field. The output probe field amplitude is 
projectively measured by a detector. 
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troduce decoherence, which has already been discussed 
extensively in the literature. In addition, we only con- 
sider one probe field, and can be easily generalized to 
multiple probe fields. 

Conditional dynamics. — At each moment, the output 
probe field 6out(i) is projectively measured by a de- 
tector, e.g., homodyne detection if the probe field is 
an optical field. We assume that (phase) quadrature 
62 = [bout{t) — 6o„t(0]/*V^ is measured with the result 
at time t being y{t). Given the measurement result, the 
plant-bath system is projected into a conditional state, 
with joint wave function lip) at t + dt given by 



m+dt)) 
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{y{t)\U{dt)\0)®\iit)). 



Here U{dt) = p-^^dt/h -g ^^^^ evolution operator; we 
assume that the input probe field (before interaction) 
is at vacuum state |0) and is separable from the joint 
plant-bath state; \y{t)) is an eigenstate of b2{t); P{y) is 
the probability density for the measurement result and 
P{y) = TYj,b{\ij{t + dt)){'il;{t + dt)\}. By integrating over 
the probe field variable, we can obtain the following non- 
linear Markovian SSE for the plant-bath state: 
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^)ak' - (flfe - al){ak' - al,)/A \tp)dt 



(2) 



and y{t)dt — —ij^k^/lki'^k — al)dt + dW/y^, which 
is from the obtained measurement result distribution 
P{y) = (dt/^)'/'exp[-(2/ + iZk^idk - with 
dW being the Wiener increment {dW^ = dt), and (ofc) = 
(tpldkltp)- When the bath is a single cavity mode, it gives 
the well-known Markovian SSE for the conditional evolu- 
tion, also known as quantum trajectory , of the plant 
and cavity mode under homodyne detection [12 1 . 

Elimination of bath modes. — To non-adiabatically 
eliminate bath modes, we apply the method by Strunz 
et al. and choose unnormalizcd coherent-state repre- 
sentation \a) = cxp[—J2k^kdl.]\0) for the bath modes. 
One can obtain an equation for \iIj{q.*)) = {a\ijj) by using 



datlipia*)) and {a[ 



4|^(a*)). The 



reduced density matrix for the plant is given by 



d a( 



■|V'(a*))(^(«)|. (3) 



gi])dt 



dpp = - -^[Hp, Pp]dt~^gk{[L\ Qk] - [L, 

k 

+ J2 VWkiQk + qI - Trp{^fc + qI}Pp)^W, (4) 



and y{t)dt = E/c\/7fc Trp[^ + Q^]dt + dW^/\/2, where we 
have introduced: 



I / d2ac-l°l'9„.|V(a*))(7^(a)|. (5) 



Qk 



Here the non-Markovianity only arises when we eliminate 
the bath, which has a memory about the plant. Eqs. (|4]) 
and ([5]) will be self-contained SMEs governing the plant 
and measurement data, if Qk can be written in terms of 
Pp and other plant operators. To derive Qk, we use the 
approach in Rcf. Q by introducing the plant operator Ok 
as follows: 



dMoc*)) = -idk{t, 



(6) 



In the simplest case, Ok does not depend on a* and 
Qk = Ok(t)pp. In general, Qk is a super-operator of pp: 

Qk = Aok(t)pp + PpAikit) + A2k{t)ppA3k{t), (7) 

where At are plant operators determined from Ok. Sys- 
tematic procedures for deriving Ok (without measure- 
ment) has been developed by Yu et al. [l3|, and applied to 
systems with different plant Hamiltonians. Yu's method 
can be generalized to our case by using interaction- 
picture \ipia*))i = U-'^{t)\4'{a*)) with a non-unitary 
evolution operator: U{t) ~ cxp[—{i/h){Hp + — 
^^^kk' ^lklk'd\.dk')t]. We can then determine Ok from 
the following consistency condition [^, [l^ : 



dt 



[dM&*))i 



da 



dt 



(8) 



In general, Ok is difficult to solve for analytically and 
must be considered case by case. In the following, we 
shall consider three interesting examples that are closely 
related to current experiments, and analytical forms of 
Ok, or equivalently Qk, can be obtained. 

Atom-cavity interaction. — As shown schematically in 
Fig. [21 we consider the following Hamiltonian: 



H =h-^az + huJcO^d + hg[(j^a) + &^a)+ 



h^[d^^{t)e 



iojQt 



(9) 



The first three terms describe the Jaynes-Cummings-type 
interaction with ujq the atom transition frequency and ct^ 



cavity 



By using the fact that / d^aafcC '"'^ |V'(q:*)) (■(/'(a)l = p. - 

/ d2ae-l°l'9a.|V(a*))(V'(a)l and from Eq.®, we ob- J- ~[ 

tain the non-Markovian SME for the plant: 
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FIG. 2: Schematics showing tlie atom-cavity system. A two- 
level atom (or a qubit) interacts with a cavity mode that 
is coupled to an external continuous optical field which is 
measured via homodyne detection. 
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the Pauli matrix, and Wc and ujq are the cavity resonant 
frequency and the laser frequency, respectively. In the ro- 
tating frame at the laser frequency, the Hamiltonian can 
be rewritten as: H = h{ujq/2)az + fi. Aa'^'a + hg{a^a^ + 
d+a) + h^J^[ah\j(t) + a^ hm{t)] with A = lOc — loq. In com- 
parison with the general Hamiltonian in Eq. ([T]) , this cor- 
responds to the case oi L = and gu — g Sik (the bath 
has only one cavity mode and we will ignore subscript k 
afterwards). By using the consisteny condition (|8]), the 
operator O = f{t)a_ and g has the following simple form: 



9 = f{t)^-P- 



(10) 



Here the time-dependent function f{t) satisfies a Riccati 
equation, / — i{LL>q — A + ij)f — gf^ — g with the initial 
condition /(O) = 0, from the assumption that the cavity 
mode is initially at a vacuum state. The corresponding 
SME for the atom density matrix reads: 



dp: 



-O-^ + .g3{/}(T+(7_, p 



dt 



— g^{f} [o"+(7-/0 + (OCT+CT- — 2 (7_(0(T+] dt + 

V^[f^-P + rp^+ - if^- + .r^+)p]dw. (11) 

This equation fully describes non-Markovian dynamics 
of the atom under continuous measurement. We can also 
obtain the corresponding master equation if we ignore 
the measurement result by averaging over dW (mean of 
dW vanishes), namely, 



(12) 



This gives the exact non-Markovian master equation for 
a two-level atom coupled to a damped cavity mode — 
a dissipative environment. Note that it differs from the 
one obtained by assuming a prior spectrum for the cavity 
mode, as have been done so far in the literature 14 1 . 

The result at the Markovian limit can be recovered 
by considering the case with the cavity decay rate much 
larger than the atom-cavity interaction rate and also the 
atom transition rate, namely 7 3> .g and ^ ^ LUq. The 
cavity mode memory becomes negligibly short and 



fit)\ 



Markovian limit 



= .9/7, 



(13) 



in which case Eq.s (fTTj) and (|T2l) reduce to the usual 
Markovian SME and master equation, respectively. 

To confirm that Eq. is the SME that correctly de- 
scribes the conditional dynamics of the atom, we numer- 
ically solve (i) the Markovian SSE for the joint atom- 
cavity wave function and (ii) the non-Makovian SME 
for the atom density matrix to see whether they both 
give the same conditional mean of ax, cry and ct^. The 
numerical results are shown in Fig.|3l We have chosen 
ujq = 1, A — 1 and 7 = 2, and the initial state for the 



atom and the cavity mode is [|-l-)2 



,]/V2(^\0). They 
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FIG. 3: The top panel shows numerical results of the time 
evolution of the conditional means: {(Jx), and (a^) given 
a particular realization of dW . The bottom panel shows the 
convergency of the accumulated numerical difference between 
the SSE and SME simulation results given different number 
of grid points for the cavity mode. 



indeed agree with each other nicely as shown by conver- 
gency of their difference. 

Linear optomechanical interaction. — We now consider 
another exactly solvable model — the linear optomechani- 
cal interaction between a harmonic mechanical oscillator 
and a cavity mode. The device is shown schematically in 
Fig.m which has been discussed extensivelyin the liter- 
ature recently [Ij. The Hamiltonian reads (l5- 17|: 



H =p^/(2m) + mujfj-^x'^ / 2 + huco) a + Hg xa^a' 



iuj{)t 



-at 6i„(i)e-'"°*] 



(14) 



Here x and p are the position and momentum of the 
oscillator with eigenfrequency uim- Since the cavity 
mode usually has a large steady-state amplitude due to 
coherent pumping by the laser, we can consider per- 
turbations around the steady-state amplitude and lin- 
earize the above Hamiltonian. In the rotating frame 
at the laser frequency, the linearized Hamiltonian is 
H = p^/2m + muj'^P/2 + hAa)a + hg' x{a) + a^) + 

^%/7[^ ^in(^) "t"^^ "whcrc g' = g a with a the steady- 

state amplitude of the cavity mode. With the same pro- 
cedure as the atom-cavity case, g can be obtained (again 
the bath has one mode with subscript k ignored): 



Q - 



{hpA^ +Ap)/{l-\h\^) 



(15) 
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FIG. 4: (color online) Schematics showing a typical optome- 
chanical device. The mechanical oscillator is coupled to a 
cavity mode via radiation pressure force. 
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withi = e-*(^-*'^)*[/o(t) + /,;(t)i] + /p(<)p. These func- 
tions /o, /i, /a: and fp are determined from the consistent 
condition, and satisfy coupled Riccati equations: 

/o = z7Tr{e + g^h - fiW - z /ig'/o/p , (16) 

=e^(^~^^)\g' + mujlfp) -ig\h + hUp) , (17) 

/, = - i(A - z7)/p - (/./™)e-^(^-^^)* - I hg' , (18) 

A = - z(A - i7)/i + .97pe'(^-'^'* - » ftg'/i/p • (19) 

These equations can be solved numerically. Similarly, if 
we average the SME over dW, we will obtain the corre- 
sponding non-Markovian master equation. It describes 
quantum Brownian motion of a harmonic oscillator cou- 
pled to a non-Markovian bath with dissipation, which 
has not yet been fully treated in the literature. 

Weak-coupling limit. — In the previous cases, we took 
advantage of the linear interaction. In general, when L 
is a nonlinear operator of the plant, there is no trans- 
parent route that leads to a closed-form solution of g. 
If the plant-bath coupling is weak, namely g^ < 7fc, we 
can perturbatively solve the problem by writing down a 
hierarchy of equations at different orders of gk/jk- The 
first-order result for the g is very elegant: 



fdT[e-^^^]kk'9k'L{^r)p 



(20) 



where L{—t) = g^-^HpT /hj^ ^iHpT /h -^-^^qj- fj-gg evolution. 

One interesting application of this result is to study 
the phonon-counting experiment recently considered in 
Refs. 
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The position of a mechanical oscillator is 
quadratically coupled to a cavity mode, namely i^int = 
hgX^{a + a') {X is the position operator normalized by 
the zero-point uncertainty). If cavity bandwidth 7 is less 
than the mechanical frequency Wm, only the time av- 
erage of — equivalent to phonon number — is impor- 
tant, and we expect a direct probe of mechanical energy 
quantization. In the proposed experiment by Thomp- 
son et al. [l3|, the coupling strength g is smaller than 
7 [i^l; we can therefore use Eq. ([20]). From X{—t) = 
X cosWrnT — P sinWrnT, we have 



9 = 



f 

Jo 



dr ( 



^X\~T)p^{g/j)Np 



(21) 



where is the phonon number, and wc have ignored 
terms proportional to e~'^* , as the characteristic measure- 
ment time scale is t ^ The resulting SME for the 
mechanical oscillator density matrix reads [cf. Eq. (j4])]: 

dp = - i[u;^N, p] dt - gcfiiX'^, [N, p]]dt 

+ ^/2^[{N, p}-2{N)p]dW + 0[{g/jf] (22) 

with (7eff = 3^/7- Note that this does not describe a quan- 
tum non-demolition (QND) measurement of the phonon 



number, as has been argued for above, since the term 
[i"^, [TV, p]] is not in the usual Lindblad form [N, [N, p]]. 
It will introduce two-phonon process and cause additional 
diffusion; we may therefore encounter unexpected fea- 
tures in the actual experiment. 

Conclusions. — We have reported a formalism that 
non-adiabatically eliminates bath modes in continuous 
quantum measurements and yields a self-contained non- 
Markovian SME for the conditional density matrix of the 
plant. Conceptually, this formalism is the mathematical 
embodiment of how memory-induced non-Markovianity 
arises when we focus on a subsystem of a larger, Marko- 
vian system. In practice, if the plant is indeed all we care 
about, the non-Markovian dynamics obtained here is an 
exact and the most efficient way of obtaining its evolu- 
tion, both in terms of analytical and numerical complex- 
ity. By averaging over measurement results, the resulting 
master equation describes the non-Markovian dynamics 
of the plant coupled to a bath that suffers from addi- 
tional dissipation, a scenario not yet fully explored in 
the literature. We have briefly illustrated the powerful- 
ness of this formalism using three examples, and we fully 
expect that it will find wide theoretical and experimental 
applications. 
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Note added. — During the preparation of this draft, we 
notice that a similar model is considered by Diosi [23| . 
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